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Abstract 

Clonality is a common characteristic of successful invasive species, but general principles underpinning the success of clonal 
invaders are not established. A number of mechanisms could contribute to invasion success including clones with broad 
tolerances and preferences, specialist clones and adaptation in situ. The majority of studies to date have been of plants and 
some invertebrate parthenogens, particularly aphids, and have not necessarily caught invasion at very early stages. Here we 
describe the early stages of an invasion by a Northern Hemisphere Hymenopteran model in three different land masses in 
the Southern Hemisphere. Nematus oligospilus Forster (Hymenoptera: Tenthredinidae), a sawfly feeding on willows {Salix 
spp.), was recently introduced to the Southern Hemisphere where it has become invasive and is strictly parthenogenetic. In 
this study, the number of N. oligospilus clones, their distribution in the landscape and on different willow hosts in South 
Africa, New Zealand and Australia were assessed using 25 microsatellite markers. Evidence is presented for the presence of 
two very common and widespread multilocus genotypes (MLGs) or 'superclones' dominating in the three countries. Rarer 
MLGs were closely related to the most widespread superclone; it is plausible that all N. oligospilus individuals were derived 
from a single clone. A few initial introductions to Australia and New Zealand seemed to have occurred. Our results point 
towards a separate introduction in Western Australia, potentially from South Africa. Rarer clones that were dominant locally 
putatively arose in situ, and might be locally favoured, or simply have not yet had time to spread. Data presented represent 
rare baseline data early in the invasion process for insights into the mechanisms that underlie the success of a global 
invader, and develop Nematus oligospilus as a valuable model to understand invasion genetics of clonal pests. 
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Introduction 

One of the major paradoxes of invasion biology is how species 
with very low genetic variation (and thus with expected low 
evolutionary potential) can still be successful invaders [1,2,3,4]. As 
invaders usually arrive in very low numbers, genetic botdenecks 
should reduce genetic diversity in the invasive range; for example: 
[5,6,7]. A number of hypotheses have been suggested to resolve 
this paradox. The first is that that genetic/genomic diversity is 
often not lost, but when it is, phenotypic variation may not be 
reduced due to plasticity [8] . In some cases, genetic diversity in the 
invasive range may increase owing to admixture between 
individuals from multiple sources [1,9,10]. Alternatively, it has 
been proposed that rapid evolution generates functional diversity 
and/or local adaptation [4,1 1]. 

One commonly observed feature of invasive species is the 
tendency to become clonal in the invasive range, which has been 
suggested to enable invasive organisms to thrive in the face of low 
diversity [2]. Clonal organisms are more than four times over- 
represented among pest invertebrates [12]. Clonality has been 



highlighted as an important mechanism by which invaders with 
low genetic diversity succeed, and a meta-analysis uncovered the 
pattern that a disproportionate fraction of aquatic invaders with 
low diversity were asexual [2]. In addition to the demographic 
benefits of asexuality (notably not investing in males, and that a 
single individual can found a new colony), this mode of 
reproduction is generally impervious to inbreeding depression 
and fit genotypes are not broken up by sexual recombination [2]. 
Conversely, while lack of the ability to shuffle pools of alleles into 
new genotypic combinations is on face value disadvantageous, 
ecological traits of parthenogens with large population sizes 
nonetheless can diversify quickly by mutation, and rapid evolution 
is an increasingly recognized contributor to invasiveness [13,14]. 

It is common for invasive organisms to transition to obligate or 
functional parthenogenesis; for example: [15,16,17,18]. These 
transitions result in the unit of selection switching instantaneously 
from the gene to a co-inherited suite of genomes (nuclear, 
mitochondrial, symbionts) at least when the transitions are 
irreversible [19]. This has important implications for the diversity 
of genotypes and phenotypes, and the fitness consequence of their 
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interactions with environmental conditions, and thus for control 
strategies [12]. An invasion that produces clones might be 
characterized by a large or small number of fixed clones, or some 

complex scenario of asexual/ sexual interactions and recurrent 
generation of new clones. Natural selection is generally expected to 
act rapidly among clonal lineages, particularly when the number 
of individuals within lineages is large, the number of hneages small, 
and they persist over time [20] . Pre-adapted genotypes should be 
favoured, while strong selection against suboptimal genot^'pes 
should occur soon after introduction. Genotvpes may be thought 
of as more or less 'generalist' or 'specialist' based on the breadth of 
their ecological niche and their environmental tolerances which 
can translate into their distributions in environmental space and 
time; both can lead to successful invasion [16,21,22,23]. New 
phenotypes can arise following mutations and karyotypic changes 
[4,1 1,24] and it is of central interest in invasion genetics to what 
extent this is an effective part of invasion ability of clonal 
organisms. 

A high proportion of studies on these topics have been of cychc 
parthenogens; in particular, Hemiptera, and more specifically 
aphids, are extremely well-represented in pest species and studies 
of them [12,25]. This important area of study requires strong 
models encompassing diverse organisms, with longitudinal studies 
starting early in invasions, ideally in multiple locations, applying 
highly-resolving genetic assays. A promising emerging model for 
understanding invasions of clonal organisms with low diversity is 
the willow sawfly Nematus oligospilus Forster (Hymenoptera: 
Tenthredinidae). This taxon has recently invaded the Southern 
Hemisphere [26]. It is native to Eurasia (from Ireland to the 
Himalayas) [27] and North America (from Alaska to Mexico) [28]. 
However, recent research points towards jV. oligospilus representing 
a group of species and the native range of the invasive species is 
unknown [26,29]. It was first detected in South America in 1980 
[30,31]. By 1993 it had spread to southern Africa [32], it reached 
New Zealand by 1997 [33] and was first detected in Australia 
(Canberra, Australian Capital Territory) in 2004 [34]. By 2009, it 
had established in most of southeast Australia (from South 
Australia to the New South Wales and Queensland border), 
Tasmania and southwestern Western Australia and was still 
undergoing expansion [35]. 

A recent study developed a high-resolution microsateUite assay 
and identified several multi-locus genotypes (MLGs) in jV. oligospilus 
in its invasive range in the Southern Hemisphere [26]. In this 
paper, we seek to advance the understanding of the evolution of 
invasive parthenogens that diverge without sexual reproduction 
and genetic recombination. In addition, we identify „'V. oligospilus as 
an exciting model taxon of parthenogenetic invasion. Here we 
assess the distribution, genetic relatedness and host associations of 
N. oligospilus MLGs in South Africa, New Zealand and Australia 
from widespread, substantial sampling ver)' early in the invasion 
process. Specifically, we address the foUowing riuestions: 1/What 
clones of JV. oligospilus are detectable in the imadcd range? 2/What 
are their geographic and host-based distributions? 3/How closely 
related are they, and can their mode of origin be inferred? 4/ Are 
there candidate 'superclones' and newly-evolved variants that 
would be valuable to monitor in time and space? 

Methods 

Field collection 

In southeastern Australia, collections occurred in December 
2007 and January 2008. The southern part of Western Australia 
was sampled in February 2010, while New Zealand was surveyed 
in February 2009. Individuals from South Africa were collected at 



one site in 2008 (n = 26). More than 1100 individuals were 
collected and genotyped from 2 1 sites in New Zealand and 69 sites 
in Australia (Table SI). All sites were on public land and sampling 
did not require permits. No endangered or protected species were 
sampled as part of this study. 

Sites were selected on the basis of presence of willow taxa, and 
distance from other sites (usually 50 to 100 km apart). At each site, 
a single mature tree per willow taxon was sampled. Foliage was 
inspected visually for sawfly presence, and any sawfly stages were 
collected into 96% ethanol, to a maximum of 30 individuals per 
tree per site. Adults were excluded because of their high mobility, 
except for South Africa due to small sample size. Defoliation of 
willow trees (i.e. when foliage on any branch had been completely 
consumed irrespectively of the amount of defoliation) was noted. 
Willow species were identified using keys for southeastern 
Australia [36,37,38] and New Zealand (van Kraayenoord et al. 
1995). 

DNA extraction and microsateUite analysis 

The DNA of up to 10 jV. oligospilus individuals was extracted per 
site per willow taxon using a salting-out DNA extraction protocol 
[39]. The type and amount of tissue used depended on the stage of 
the sawfly. For adults, DNA was extracted from the abdomen, 
while for larvae and pupae, the head and a portion of the thorax 
were used. Due to their small size, whole first-instar individuals 
were used. 

Twenty-five microsateUite loci were used to identify multilocus 
genotypes (MLGs) that differ by one or more allele [26] . Fourteen 
of the loci were polymorphic in the Australian, New Zealand and 
South African samples, ten of which contributed to discrimination 
among genotypes. The loci that did not discriminate clones were 
nonetheless considered worth including (particularly given the 
convenience of simultaneous amplification of many loci in this 
system) for potential detection of rare/ new variants since they are 
polymorphic in ,\. oligospilus in locations with males as well as 
females, and in other tenthredinid taxa [26]. Use of 25 loci should 
yield high resolution: previous studies have attributed clonal 
identity in introduced populations on the basis of many fewer loci, 
e.g. 5 to 7 microsateUite markers [40,41]. Individuals bearing the 
same 25-locus genotype wiU be referred to as members of the same 
clone, whUe acknowledging that even mother-daughter pairs are 
likely to differ by at least a few point mutations somewhere in the 
genome. One primer of each primer pair was labeUed with infra- 
red dyes IRD700 or IRD800 (LI-COR Biosciences). Using die 
Qiagen PGR multiplexing kit following the manufacturer's 
instructions, 12 and 13 loci were multiplexed in 10 |J,1 reactions, 
with a maximum of seven loci per channel (Multiplex A: 1RD700: 
WF2, WFll, WF15, WF21, WF26, WF27, WF40; IRD800: 
WF12, WF13, WF25, WF30, WF39, WF41. Multiplex B: 
IRD700: WF9, WF20, WF23, WF32, WF34, WF35, WF38, 
IRD800: WF6, WFIO, WF14, WF16, WF37). PGR products were 
electrophoresed on 6% acrylamide gels in LI-GOR Global 1R2 
two-dye DNA sequencers (models 4200 and 4300). MicrosateUites 
were viewed in LI-COR software Saga Generation 2.0, and scored 
visually. 

Analyses 

Genetic analyses we're conducted according to sampling units 
based on country (Australia, New Zealand and South Africa), and 
within countries, regions separated by major land or water barriers 
between populations of jV. oligospilus. The rationale for the choice 
of these geopolitical units is that their physical organization, 
economic activities and biosecurity implementation are likely to be 
dominant factors in the likelihood of introduction. Six regions 
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were defined: two in New Zealand (the North and South Islands), 
three in Australia (southeastern Australia; southwestern Western 
Australia isolated from southeastern Australia by ~1 100 km of the 
NuUarbor Plain, and Tasmania isolated by ~400 km from the 
nearest part of mainland Australia by the Bass Strait), and South 
Africa. Genotypic diversity (G/N) was calculated for each country 
and region, using the number of genotypes (G) divided by the 
number of individuals genotyped (N). Shannon-Weaver diversity 
index, allelic diversity, observed (Ho) and expected heterozygosity 
(He) were calculated in GenAlex 6.4 [42]. Shannon- Weaver 
diversity index [43] was expressed as to allow comparison based 
on the number of genotypes present in the populations. 
Departures from Hardy-Weinberg equilibrium (HWE), linkage 
disequilibrium and f-statistics were calculated and assessed with 
exact tests in Genepop 4.0 [44]. Pair-wise Fst estimates between 
populations were calculated based on the 25 loci with combined 
statistical significance calculated according to Fisher's method as 
implemented in Genepop. HWE analyses were performed and 
expected and observed heterozygosity calculated initially for the 
full dataset, and then for a reduced dataset containing one 
individual per genotype to avoid distortions caused by multiple 
representatives of the same clonal genotype [40,45]. 

Unrooted neighbour-joining trees were constructed 
based on Gavalli-Sforza's chord measure of genetic 
distance [46] in Populations 1.2.32 (www.bioinformatics. 
org/~tryphon/populations/#ancre_fonctionnalites). One tree 
was built to depict the relationships among MLGs (inferred allele 
changes were mapped onto branches), and one tree depicted the 
population network on the basis of one individual per MLG per 
region. To assess genetic variance within and between countries 
and regions, a nested Analysis of molecular variance (AMOVA) 
was carried out with regions nested within countries in GenAlex 
6.4 [42]. P-values were assessed with 999 random permutations. 

Results 

Detection and distribution of clones 

The combination of 25 microsateUite markers discerned 16 
different MLGs in the Australasian and South African invaded 

range sampled (Table 1). Of the MLGs identified, ten were found 
in Australia, seven of which were unique to that continent. 
Similarly, seven MLGs were found in New Zealand, five of which 
were unique, while in South Africa, one of three clones was unique 
(Table 1, Fig. 1). The 'Main clone', the most widespread MLG, 
occurred over most sampled areas of Australia and New Zealand 
and was also present in South Africa (Figs. 1, 2, 3). It represented 
65.7'l'o of all individuals genotyped. It was, however, absent from 
Western Australia. The second most widespread MLG, the 'North 
clone' comprised 24.9% of individuals. It encompassed the 
northern part of southeastern Australia as well as most of Western 
Australia, although it was present in low proportion at several 
other locations (e.g. one site in Tasmania) (Fig. 3). It was also 
found on North and South Islands of New Zealand. It appeared to 
be spreading, because it was found in high proportions at some 
sites but smaller numbers in surrounding sites (Fig. 2). Although 
both clones were widespread, they rarely overlapped within a site. 

Other than 'Main' and 'North' clones, most MLGs were 
localized and had low frequencies. Three exceptions were clone 
Aus 1 , which was found in several sites in Western Australia around 
the port city of Albany and was unique to this region, and Aus2 
and AusS which were common locally (Table 1, Fig. 2 and 3). 
Clone Aus2 occurred at one site, along with the Main clone. Clone 
AusS occurred at two isolated sites, where it was the only clone 
present (Table 1, Fig. 3). Other MLGs were usually represented by 



one individual per site. While two MLGs were found in more than 
one site (Aus2, and Aus 4), other MLGs were unique to a single 
site. Locally common MLGs were found on all willow taxa within 
a site. Few statements about distributions can be made for South 
Africa since all individuals were sampled at one site. However, one 
rare MLG was present in South Africa and South Australia (Aus6). 

The two most widespread MLGs, 'Main' and 'North' clones 
were found on all willow taxa sampled showing no obvious host 
specialization. Similarly, locally abundant MLGs Ausl and Aus2 
were found on all the willow taxa present at the sites where they 
occurred. Only the 'Main' 'North' and Ausl MLGs have been 
found at sites where willows suffered defohation. 

Clone relatedness and genotypic diversity 

The 16 MLGs sampled were extremely similar at their 25-locus 
genotypes: they differed by a maximum of 3 alleles and 3 loci 
(Table 1). Although other interpretations are possible, the data are 
most parsimoniously consistent with most or all of the 15 MLGs 
being derived by mutation from the most common and 
widespread Main MLG. On an unrooted neighbour-joining tree 
of the MLGs, 'Main' is the most interior, and from it branch eight 
individual MLGs differing from Main each by a single allele 
difference (Fig. 4). Another three branches from Main each 
contain two or three MLGs linked by progressive single 
differences. One of the multi-MLG branches was limited in 
distribution to New Zealand (NZl and NZ2). The other two 
branches were not geographically restricted; one encompassed two 
locahzed MLGs found in Australia (Aus5 from Western Australia, 
Aus7 from Southeastern Australia) and also 'North', found in 
Australia and New Zealand, while the last multi-MLG branch 
spans Ausl from W(;stcrn Australia, and SA from South Africa. 
While the tree-building process did not allow reticulation, none is 
necessary to explain the data: MLG relationships can be depicted 
by simple sequential changes without homoplasy or ambiguity. 

Overall, the three countries and regions had low genetic 
diversity based on genotypic diversity, (P and allelic di\'ersity 
(Table 2). Shannon-Weaver diversity index was extremely similar 
for the three countries and regions, between 1.35 and 1.37. 
Observed heterozygosity and allehc diversity were equally low in 
the regions, respectively ranging from 0.42 to 0.44 and 1.48 to 
1.68. There were frerjuent deviations from HWE in the dirertion 
of heterozygote excess, a common property of microsateUite data 
sets from obligate parthenogens proposed to diversify by mutation 
[47] (Table 2). The number of private alleles was low: New 
Zealand had three private alleles, while Australia had four and 
South Africa had none (the latter result likely contributed to by 
small sample size). 

Low genetic differentiation among countries and regions 

-FstS for all country or region pairs were low, ranging from less 
than 0.001 to 0.05. All pairwise comparisons showed that sampling 
locations (countries or regions) did not differ significantly in allele 
frequencies (p>0.05). According to the AMOVA, most of the 
molecular variance, 88%, was attributed to within-region effect, 
while only 12% was attributed to between-region effects and less 
than 1 % to between-countr)' effects (Table 3). The population tree 
placed the regions into two groups. The main grouping included 
southeastern Australia, Tasmania, the North and South Islands of 
New Zealand, with South Africa very slightiy outside that group, 
and with Western Australia placed on a long branch slightiy more 
allied to South Africa (Fig. 5), presumably owing to the close 
relationship of Western Australian-limited MLG Aus 1 and South 
Africa-limited SA. 



PLOS ONE I www.plosone.org 



3 



June 2014 I Volume 9 | Issue 6 | e97744 



Superclones May Explain Invasion Success 



g g 



ggggggggg 



e c 

O is 



:2 £ :2 



S > 



Dl Ol Oi 



a S a 



^ 



a 



o 

01 



< 

g 



< 

g 



5 

B 2 

c .™ 

3 = 



£ 

< 



° ° E c o 

u u o 

£ ■£ £ S E 



o 



N N NJ M N 
? ? f J ? 



19 
-I 



N N N M tsl 



PLOS ONE I www.plosone.org 



4 



June 2014 | Volume 9 | Issue 6 | e97744 



Superclones May Explain Invasion Success 





PLOS ONE I www.plosone.org 



5 



June 2014 | Volume 9 | Issue 6 | e97744 



Superclones May Explain Invasion Success 



a) 



b) 




MLGs 

□ Main 

■ North 

■ NZl 

□ NZ2 

■ NZ3 

□ NZ4 

□ NZ5 

■ AUSl 

□ AUS2 

■ AUSB 

■ AUS4 

□ AUS5 

■Ause 

■ AUS7 

■ AUS8 

■ SA 



Figure 3. Frequency of multilocus genotypes (MLGs) of Nematus oligospilus in Australia. Each circle represents a location sampled and the 
relative frequency of each MLG within the site. Each color represents a different IVILG. a) southeastern Australia and Tasmania b) Western Australia. 
doi:10.1371/journal.pone.0097744.g003 



Discussion 

Clone distribution and relatedness 

Nematus oligospilus is a very successful parthenogenetic invasive 
organism with low genotypic variation, and two extremely 
successful genotypes. The 'Main' MLG has colonized three 
continents and was present in all regions studied except Western 
Australia, while the 'North' MLG was widespread in Australia, 
and New Zealand. They were both found on all willow taxa 
present, indicating little host specialization. This, coupled with 
their wide geographic distribution and the very high proportion of 
individuals sampled belonging to them, means that they fit 
published definitions of 'superclones': genotypes dominating 
populations over large areas; for example: [41,48]. Other clonal 
invasive organisms also have only a few widespread genotypes in 
their introduced range. Examples are relatively common in plants. 



with many cases of one or a few dominant clones in their invasive 
populations; for example: [22,49] and in some invertebrates such 
as aphids [21,40,41,48]. Our data are unusual in that they are 
particularly detailed (25 loci compared to as few as five 
microsatellites in similar studies [41]) and samples animals from 
early in the invasion process. The widespread distribution of the 
two very common MLGs, their wide host range and the fact they 
were found singly at sites where willows had been defoliated may 
indicate that characteristics of these clones contribute to early 
invasion success. Our study expands the published literature to 
encompass highly-resolving genotypic patterns in a representative 
of the family Tenthredinidae, a dominant component of the 
Hymenopteran parthenogenetic fauna [12]. 

Clones were closely related, suggesting that the number of 
introductions into the Southern Hemisphere has been few, except 
perhaps if they came from a stock of unusually limited diversity. 
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NZ2 




Aus3 



Figure 4. Unrooted neighbour-joining tree based on Cavaili-Sforza cliord distance depicting relationships between different 
multilocus genotypes of Nematus oligospilus. Under the assumption that multilocus genotypes (MLGs) more interior in the networic are the 
ancestors of more exterior ones on the same branch, inferred mutational changes are depicted on the branches, e.g. Wf21 (337^333) indicates that 
MLG NZ5 could have arisen most parsimoniously by a 4 base pair deletion of the Wf21^^^ allele from an ancestor with the genotype of Main MLG. 
dor:1 0.1 371/journal.pone.0097744.g004 



We infer on the basis of the network of MLGs, with support from 
geographic and numerical distributions, that many of the clones 
are descended by mutation from an early invader (very likely the 
'Main' clone). For example, some of the rarer clones seem to have 
arisen in situ, as they occurred in very low proportions at sites 
where a common MLG predominated and differed from it at only 
a single allele at a single locus. Following introduction in new 
environments, clonal lineage changes can occur through mutations 
[4] as observed in other studies of invertebrates and plants 
[11,49,50]. The dominance of localized MLGs could be due to 
being favoured by environmental conditions. Conclusive determi- 
nation of whether these MLGs were introduced separately or have 
arisen in situ would require further work, but given their 
geographic abundance distributions and heterozygous excess (a 
characteristic of mutational divergence in clones, [47]), we 
currently favour the latter explanation. Possibly, rarer MLGs 
have been or will be favoured by environmental conditions; such a 
process could explain the presence of some localized MLGs 
dominating some site. However, the present paper presents one of 
the few cases globally where such a process could be inferred, and 
so N. oligospilus and the newly-available genetic tools for it offer 
many possibilities to test specific hypotheses concerning the success 
of clonal invaders [51,52]. 

Host association 

Host plants can provide selection pressures that lead to 
speciahzation of herbivore populations [53,54]. While clonal 
specialization has been observed in some invading aphids [16,45], 
the widespread MLGs of M. oligospilus did not depend on 
association with specific willow taxa, indicating that this species 
may be a host generalist. However, in the field, when several 



willow hosts are present, population densities of jV. oligospilus are 
constandy higher on some host taxa than others [35]. The less 
favoured hosts may be more challenging and exert higher selection 
on N. oligospilus and increase the chance for mutations to spread. In 
a study on Sitobion aphids, although no host specialization was 
found, rarer clones had enhanced performance on the most 
chemically defended host plants, and it was hypothesized that 
rarer clones were relatively favoured partly because of the 
presence of chemical defences [55]. Due to the low numbers of 
the rarer clones, this could not be assessed in this study, but 
warrants future investigation. 

Introduction and dispersal in the landscape 

The native range source of the N. oligospilus that have invaded 
the Southern Hemisphere has not yet been resolved, and genetic 
data can contribute to clarification of relationships of these insects 
around the world [29,35]. The modest number of extremely 
closely-related clones, all derivable from the 'Main' clone by very 
few mutations, indicates a small number of introductions from a 
single source. AU Australasian and the Kmited number of South 
African individuals assayed are possibly derived from a single 
founding clone. However within the sampled invaded range, 
multiple introductions may have occurred. For example, the 
unique clone in Western Australia ('Aus 1 ') is more closely allied to 
a South African MLG than it is to Australian clones (although the 
relationships are all extremely close). Thus Ausl may represent an 
independent introduction into Australia, quite possibly from South 
Africa, although a more thorough sampling of the African 
continent is required. Many unwanted introductions to Western 
Australia have occurred from South Africa, including many plants 
and invertebrates such as the African black beede, now considered 
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South Island (NZ) 




Western Australia (A) 



North Island (NZ) 
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Figure 5. Unrooted neighbour-joining tree based on Cavalli- 
Sforza cliord distance depicting relationships between differ- 
ent regions where Nematus oligospilus occurs in the Southern 
Hemisphere (SE: Southeastern, A: Australia, NZ: New Zealand) 
on the basis of the multilocus genotypes (MLGs) present in 
each region. 

doi:10.1371/journal.pone.0097744.g005 
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an important pest [56,57]. An independent introduction into 
Western Australia is supported by reports of sawflies first being 
seen there five years later than in the rest of Australia (2004) (F. 
Ede, unpublished). 

Sawflies in general are not strong flyers [58]. Therefore, self- 
dispersal could occur over limited areas but is unlikely over large 
water- and land-barriers. For example, it is highly unlikely that N. 
oligospilus flew from South Africa to Australasia, although we 
cannot rule out transport by wind. Similarly, it seems improbable 
that 'North' MLG found in southeastern Australia and Western 
Australia dispersed unassisted across 1 1 00 km of the NuUarbor 
Plain (arid habitat lacking willow hosts), or that the single 
individual of the 'North' clone sampled in Tasmania self-dispersed 
to the location of its capture close to Launceston airport. Most of 
the dispersal seen between regions is more likely due to human- 
mediated jump dispersal, which is also a major cause of long- 
dispersal in many other invasive species; for example: [59,60]. 
When ready to pupate, larvae move or drop to the ground, 
spinning cocoons on vegetation or in soil litter. In high infestation 
areas, cocoons can be found on any surface available including 
outdoor furniture and cars (Caron, pers. obs.). Cocoons could 
therefore easily be transported inadvertently for long distances. 
Furthermore, in at least a few cases, people have moved sawflies 
deliberately in the hope of controlling willow infestations (F. Ede, 
pers. obs.). Human-related dispersal of invasive organisms is 
common and observed in other introduced sawflies [61] and other 
parthenogenetic pests such as grape phylloxera [Daktulosphaira 
vitifoliae Fitch), in which dispersal among localities is due solely to 
human-mediated movements of infected plants [62]. However, 
self-dispersal plausibly explains much of the spread of the Main 



PLOS ONE I www.plosone.org 



8 



June 2014 | Volume 9 | Issue 6 | e97744 



Superclones May Explain Invasion Success 



Table 3. Results from the nested Analysis of Molecular Variance (regions nested within countries) (df: degree of freedom, SS: sum 
of squares, MS mean squares. Est. Variance: Estimate of variance, % variation: percentage of total variation. 





Source 


df 


SS 


MS 


Est. variance 


% variation 


P 


Among countries 


2 


66.006 


33.003 


<0.001 


0% 


1.0 


Among regions 


3 


128.705 


42.902 


0.325 


12% 


0.001 


Within regions 


1151 


2770.058 


2.407 


2.407 


88% 


0.001 



P: P-value (calculated from 999 permutations). 

doi:1 0.1 371 /joumal.pone.0097744.t003 



and North MLG at smaller spatial scales within southeastern 
Australia and New Zealand. 

Conclusions 

Two dominant MLGs (or superclones) of jV. oligospilus occur in 
the Southern Hemisphere invasive range. Empirically these 
superclones are highly successful invaders, comprising most of 
the spread of jV. oligospilus through the colonized range. Whether 
this success persists over longer periods remains to be tested. 
Clones with broader preferences and tolerances may be better at 
colonizing new areas, but may not be as good at competing [63]. 
Thus we expect invasions to be dynamic in their demogenetics: for 
example, one or a few genetic variants can be closely associated 
with dispersal abihty but antagonistic to reproductive success [64] . 
For jV. oligospilus, new MLGs seem to have arisen in situ within the 
Southern Hemisphere. Although two main MLGs now dominate, 
new and rare MLGs may be favoured subsequently, following 
changes in environmental conditions and as non-equilibria such as 
founder effects dissipate [3] . Furthermore, additional introductions 
of distinctive MLGs to the introduced population could provide 
new impetus to invasion, as seen with the multiple introductions of 
specialised clones of the pea aphid [16]. While few introductions 
seem to have occurred in the present case, the extreme success of 
N. oligospilus in the Southern Hemisphere highlights the impor- 
tance of biosecurity. The introduction of a single individual 
capable of reproducing asexually and pre-adapted to the 
conditions of the new environment could lead to a successful 
invasion. 



We propose that Nematus oligospilus is a valuable model in 
invasion genetics. Data and tools have been established early 
during invasion of significant sections of the globe. We encourage 
repeated and broad sampling of the invaded range and potential 
native range to help infer the mode and location of generation of 
MLGs. Coupled with detailed studies of the attributes of tractable 
individual clones, this can gready add to the understanding of 
invasion genetics of clonal pests. 

Supporting Information 

Table SI Sites and willow taxa surveyed. N: number oi Nematus 

oligospilus genotyped. 

(DOCX) 
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